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ABSTRACT 

Complex organic molecules are ubiquitous companions of young low-mass protostars. Recent 
observations suggest that their emission stems, not only from the traditional hot corino, but 
also from offset positions. In this work, 2D physicochemical modelling of an envelope-cavity 
system is carried out. Wavelength-dependent radiative transfer calculations are performed and 
a comprehensive gas-grain chemical network is used to simulate the physical and chemical 
structure. The morphology of the system delineates three distinct regions; the cavity wall 
layer with time-dependent and species-variant enhancements; a torus rich in complex organic 
ices, but not reflected in gas-phase abundances; and the remaining outer envelope abundant 
in simpler solid and gaseous molecules. Strongly irradiated regions, such as the cavity wall 
layer, are subject to frequent photodissociation in the solid phase. Subsequent recombination 
of the photoproducts leads to frequent reactive desorption, causing gas-phase enhancements 
of several orders of magnitude. This mechanism remains to be quantified with laboratory 
experiments. Direct photodesorption is found to be relatively inefficient. If radicals are not 
produced directly in the icy mantle, the formation of complex organics is impeded. For ef¬ 
ficiency, a sufficient number of FUV photons needs to penetrate the envelope; and elevated 
cool dust temperatures need to enable grain-surface radical mobility. As a result, a high stellar 
luminosity and a sufficiently wide cavity favor chemical complexity. Furthermore within this 
paradigm, complex organics are demonstrated to have unique lifetimes and be grouped into 
early (formaldehyde, ketene, methanol, formic acid, methyl formate, acetic acid, glycolalde- 
hyde) and late (acetaldehyde, dimethyl ether, ethanol) species. 
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1 INTRODUCTION 


Young, forming stars emerge in regions seeded with diverse 
molecules . From highly abundan t deuterated species in prestellar 
cores tsee ICeccarelli et al.ll2014l for a revi ew) to complex organic 
molec ules in hot cores and corinos (see iHerbst & van DishoeckI 
I 2 OO 9 I for a review), chemical variety spans the entire lifetime of 
a protostar. Complex organic species are particularly alluring due 
to their potential astrobiological implications. 

Complex organic compounds are loosely defined in both 
chemistry and astronomy, b ut typically mean large Q 6 atoms) 
carbon-containing species (iHerbst & van PishoecM I2OO9I) . They 
were first o bserved in hot cores surroundi n g high-mass pro- 
tostars (e.g., iLovas. Johnson & Snvdej 1 19791 : Blake et ^ 1 19871 : 
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[Kuan & Chuang||2014h . but have since also bee n detected in the en- 


virons of several low-mass counterparts (e.g.. van Dishoeck et al 

1995 

Cazaux et alj|200 

3: Bottinelli et alJl2004allbl: J0reensen et al 

2012 

Sakai et alJl201Sl; 

Maurv et alj|20l4). The oroximitv of low- 


mass protostars facilitates spatially resolved studies. So far, single¬ 
dish telescopes have not revealed the origin of complex organics 
conclusively, but this is changing - progress is being made with 
powerful interferometers, such as the Submillimeter Array (SMA), 
the Plateau de Bure Interferometer (PdBI) and the Atacama Large 
Millimeter/submillimeter Array (ALMA). 

It is well established that hot corinos (defined as the in¬ 
ner warm zones where Tdust ^ 100 K) are chemically rich; 
however, several studies have shown that colder envelopes and 
positions with impinging outflows may also glow with emis¬ 
sion from complex organic molecules. Formic acid (HCOOH), 
methyl formate (HCOOCH 3 ), acetaldehyde (CH 3 CHO), methyl 
cyanide (CH 3 CN) and ethanol (C 2 H 5 OH) were detected at a dis- 
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tance of ~ 10"* AU from the L 1157-mm low-mass protostar 
11 Lq, ISugimura et alj I 201 II) - towards the B 1 position in 


the brightest carbon monoxide (CO 

clum 

p in the blue outflow 

lobe (Arce et al. 20081: ICodella et alJ 

0 

0 

Sugimura et alj|201ll; 


ready present on the grain surfaces and are subsequently liberated 
as gases by passing shocks, while others (e.g., methyl cyanide and 
acetaldehyde) form in the gas phase. 

An outflow-affected position in the Perseus Bl-b dust core 
(near a namefellow ~ 3 Lq protostar) and another in the Ser¬ 
pens core (specifically, SMM4-W, separated by ~ 3000 AU 
from SMM4, which has a luminosity of~ 5 Lq) were studied by 
lOberg et"^ ( l2010ah and lOberg et alj i201lh . respectively. The for¬ 
mer did not show emission from complex organic species in con¬ 
trast to the latter, which was bright in acetaldehyde and dimethyl 
ether (CH 3 OCH 3 ). The envelopes of Bl-b, SMM4, and SMMl 
(~ 30 Lq) contain lines from ketene (H 2 CCO), formic acid, 
methyl formate, acetaldehyde, dimethyl ether, methyl mercaptan 
(CH 3 SH), propynal (HCCCHO), r nethoxy (CH 3 C )) and potentially 
oxirane (CH 2 OCH 2 ) iOberg et alj|2010a . 1201 ll : ICernicharo et alj 
l2012h . The envelope encompassing the most luminous target, 
SMMl, has the highest overall abundance of complex molecules. It 
was proposed that the ice content was revealed thanks to photodes¬ 
orption, which may be efficient on larger scales due to the presence 
of o utflow cavities, thr ough which FUV photons can escape. 

Ijaber et al] (1201 4h presented data from The IRAS 16293 Mil¬ 
limeter and Submillimeter Spectral Survey (TIMASS. ICaux et alj 
I 2 OI ih . Detections of ketene, methyl formate, acetaldehyde, 
dimethyl ether, methyl cyanide and formamide (NH 2 CHO) in the 
cold (< 50 K) envelope (spanning ~ 6000 AU) surrounding 
IRAS 16293-2422 A and B (combined 22 Lq) were asserted. 
This is supported by earlier work on this source, which showed spa¬ 
tial segregation for various species (e.g..l\^ Dishoeck et alJll995l: 


lCecc^elh^tayj2QQ^ Schbier et al.ll2()02l : Bisschop et alj|2008h . 


Oberg et ^ j201lh. based on t heir observations and labora¬ 
tory experiments dOberg et alj|2009h . argued that complex organ¬ 
ics form sequentially. While CO is frozen out, HCO-rich chemistry 
prevails, leading to ‘cold’ species like methyl formate and acetalde¬ 
hyde. Upon warm up, CH 3 / 2 -rich molecules are formed, building 
‘hot’ species like dimethyl ether. However, observations of ketene, 
methyl format e, acetaldehyde, and d imethyl ether in the prestellar 
core L1689B dBacmann et alj|2012h. and form ic acid towards the 
dark cloud L183 1 Reouena-Torres et alj|200'^ . argue against this 
sequential chemistry scenario due to the lack of elevated tempera¬ 
tures in these cold (~ 11 — 12 K) sources. 

A third prestellar core, L1544, was observed to show emission 
from ketene, formic acid, acetaldehyde, and propyne (CH 3 CCH), 
but not from other species (s uch as methyl formate and dimethyl 
ether) found towards L1689B jVastel et alj2014l) . ICernicharo et alj 
( 1201 2 h secured a detection of the rnethoxy radical towards Bl-b in 
Perseus and, based on laboratory data, argued that it is most likely 
produced in the gas phase from reactions between CH 3 OH and OH 
- a reaction which has since been calculated to have a significantly 
larger rate coefficien t at low temperatures than previously thought 
dShannon et alj20l^ . In summary, complex organic molecules are 
associated with outflow cavity walls, envelopes and even cold, dark 
cores, but it remains unclear when and how the various species 
form. 

The complete chemistry of complex organics has been elud¬ 
ing astrochemists. Models have shown that gas-phase pathways 
contribute to the complex organic budget only under hot core¬ 
like conditions, when the densities and gas temperatures are high 


(~ 10 ’*' cm~^ and >100 K , e.g.. lChamlev. Tielens & MillaJ 19921: 
[Rodgers & Charnie^l200lh . Instead, complexity could be built up 
in the ic es. One solid-phase scenar io is based on atom addition re¬ 
actions dTielens & Chaml^ll997ll . This has been experimentally 

J iroven for the sequentia l hydrogenation chain l eading to methano l 
Tielens & Hagenl|l982l : IWatanabe et alj[^004l : iFuchs et^l2009h . 

An alternative scenario is that of processed ‘zeroth generation’ ices 
(i.e., simple ices like CO, methane (CH 4 ) and CH 3 OH), which 
have been subjected to far-ultra violet (FUV) irradiation, heatin g 
and/or cosmic rays (CRs, e.g.. iGarrod. Weaver & HerbstI 1200^ . 
Experimentally, an array of species is produced upon F UV irra¬ 
diation of CH 3 OH and CO mixtures dOberg et al.ll2009h . Recent 
ideas incl ude combinations of radical-radical and hydrogenation 
reactions dFedoseev et alJl^lSh . and revisions of gas-phase reac¬ 
tion rates, potentially making gas-phase chemistry non-ne gligible 
at lower temperatures dBalucani. Ceccarelli & Taouetl20151) . In ad¬ 
dition, alternative mechanisms, like Eley-Rideal and complex- 
induced reactions dRuaud et al.l l2015h . or CR-induced diffusion 
dReboussin et alj|2014h . have also been invoked to pave the way 
to complexity. However, efficient mechanisms that couple ices to 
the gas phase resulting in appreciable observed abundances remain 
challenging. 

To obtain more insight, this paper offers physicochemical 
models of envelopes harboring low-m ass protostars, i ncluding the 
impact of outflow cavities. Previously. [Bruderer et al.l d200(Al201^ 
modelled a high-mass young stellar object (YSO) with an outflow 
cavity with a focus on ions and diatomic hydrides. The authors 
pre-calculated chemical abundances with a gas-phase network at 
a given time, which saves computational time. Instead, here, a two- 
phase chemical netwo rk is used for time -dependent calculations 
on a cell-bv- cell basis. [ Viss er et^ d2012h simulated the low-mass 
analogue of [Bruderer et al. d2009L |201(]|) . but focused on a thin 
warm/hot layer along the cavity wall, from which the far-infrared 
CO and water emission originate, rather than the larger cooler lay¬ 
ers of the envelope that are central to this work. 

In this work, the dust temperature and the stellar radiation held 
of the system are computed with RADMC - a 2D wavelength- 
dependent radiative transfer module. A comprehensive gas-grain 
chemical model is employed in order to simulate the distribution 
of complex organic molecules across the envelope as a function of 
ti me. The major upgr ade in this paper, in comparison to the setup 
of ivisser et al.l ( l2012li . is the inclusion of various grain-surface pro¬ 
cesses. Our models test the hypothesis that enhanced irradiation 
of the envelope, due to outflow cavities, stimulates the observed 
appearance of complex organic molecules in the gas phase in the 
colder offset positions from the protostar. The physicochemical 
models are described in Section The results are presented in 
Section [ 3 ] and discussed in Section |4] The main findings are sum¬ 
marised in Section|5] 


2 PHYSICOCHEMICAL MODEL 
2.1 Physical setup 

The goal of this work is to model the physical structure and the 
chemical composition of the envelope affected by outflow cavities 
on scales of 10 000 AU. A static setup is adopted, representative of 
a low-mass embedded protostar (Class O/I). The chemistry is then 
computed in time at each point under constant physical conditions. 

The gas density distribution (nn) sets the underlying dust dis¬ 
tribution upon the assumption of a gas to dust mass ratio (rgd). 
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Table 1. Density distribution, outflow cavity, and stellar param¬ 
eters (fiducial model setup). 


Parameter 

Units 

Value 

P 


-1.7 

Hn 

AU 

35.9 

^H 2 ; in 

cm"® 

4.9 X 10® 

ttcav 

AU 

4.3 X 10"* 

bcav 

AU 

6.5 X 10® 

a{z = lOOOAU) 

deg 

109 

a (2 = 10 OOOAU) 

deg 

45 

T* 

K 

10 000 

L, 

L© 

35.7 


Once the dust distribution is coupled with a set of stellar vari- 
ables, the temperature and radiation field di stribut ions are defined. 
Jo rgensen. Schoier & van Dishoeckj ( l2002h and iKristensen et alJ 
I2OI0I) fitted dust cont inuum emission stemm ing from envelopes 
with the Dusty code dlvezic & Elitzui|[l9^ . DUSTY models are 
spherically symmetric power law density distributions. The authors 
imposed the spectral energy distributions (SEDs) and continuum 
radial profiles at 450 and 850 )tm from the JCMT as restrictions 
on these models to arri ve at a set of bes t -fitting parameters. Those 
results were updated in iKristensen et al] ( 120121) based on new Her- 
schel fluxes. In this work, the density distribution derived with the 
Dusty model of NGC 1333-IRAS2A is taken as a test case for 
low-mass protostars with stron g emission from c omplex orga nics. 

In the modelling work of Kristensen et al.l j2O10Ll2Ol^ . out¬ 
flow cavities were neglected. Here, their presence is simulated by 
carving out ellipsoidal cavities from the DUSTY d ensity distribu¬ 
tion. This ha s pr eviously been done in the work of ISruderer et alj 
ll2009Ll2010l) and IVisser et al.l (1201 2h with two slightly different pa- 
rameterisations. In this work, equation 1 from ivisser et al.l ( 1201 2h 
is used. 


R 


cav 




( 1 ) 


where (i?cav, ^cav) are the cylindrical coordinates of the cavity wall, 
and flcav and bcav are the semimajor and semiminor axes, respec¬ 
tively. Within the outflow cavity a comparatively low constant den¬ 
sity is assumed, 

TiH (cavity) = 2 x 10'*cm~^, (2) 


and elsewhere. 


riH (elsewhere) = nu; in x 



(3) 


where r is the radial spherical coordinate, p is the power index 
of the density profile, rin is the inner boundary of the envelope 
and nu; in is the H nuclei number density at rin. It is assumed that 
riH = 2 X riHa (at all positions). With the selected values of 
and 6cav, the full opening angle of the out flow cavity is 109° a t 
2 = 1000 AU and 45° at z = 10 000 AU. IPlunkett et alj <20131) 
tabulated the full opening angles for this source (table 4, although it 
is unclear at what distance from the source the values are derived), 
and find 30° and 56° for the two outflows in the blue lobe, and 
32° and 93° for the two in the red lobe. These values are compa¬ 
rable to the full opening angles in this model. For other sources, 
wider angles have also been mea sured, e.g., 110° for HH46/47 
dVelusamv, Danger & Marshll2007l) . 

The amount of radiation stemming from a star depends on its 


blackbody temperature and the size of the emitting surface, i.e., on 
the stellar temperature (T*) and radius (Rt ). The bolometric lumi¬ 
nosity of NGC 133 3-IRAS2A is 35.7 Lq ( IKristensen et~Zll2012l: 
iKarska et alj|2013h . Protostars are thought to have an ultraviolet 
(UV) excess due to accre tion of material from the disc onto the star 
(e.g.. lSpaans et al.lll995l) . FUV photons are particularly efficient at 
photodissociating and photodesorbing molecules. Therefore in this 
work, a stellar temperature of 10 000 K is adopted in order to ac¬ 
count for the anticipated UV excess. From the definition of stellar 
luminosity (L*), the following parameterization is used to obtain 
the stellar radius: 


R, 


(T./Tq)^ 


R. 


©• 


(4) 


Although parameters derived for NGC 1333-IRAS2A are used, our 
aim is not to model this particular source, but inst ead to build a tem- 
plate for a Class 0 protostar (see, e.g., fig. 13 in [Evans et alj|20()^ 
for statistics per Class). The set of adopted parameters is summa¬ 
rized in Table[T] 

Once the dust de nsity distribution and stellar p roperties are 
prescribed, R ADMcQ dPullemond & Dominikll2004l) is employed 
to calculate the temperature and radiation field distributions. An rgd 
of 100 is assumed and the dust mass density (pdust) is given by: 


Pdust = 0.5 X nu X p X mp/rgd, (5) 


where rrip is the proton mass and p is the mean molecular mass 
of the gas (taken to be 2.3). The dust is assumed to be a mixture 
of carbonaceous material (25 per cent) and silicates (75 per cent). 
Opacity tables for bare g rains covering a range of grain sizes from 
[Pontoppidan et al.l <20071) are adopted. Modifications of the opac¬ 
ities to include icy mantles are expected to alter the dust temper¬ 
atures (Tdust) by, at most, a few K (M.K. McClure, priv. comm.). 
Computations are carried out over frequencies in the 2.998 x 10^° — 
3.288 X lO*^® Hz (0.09— 1000 pm) range, ensuring coverage of the 
dust’s crystalline silicate features and the FUV range. The FUV ra¬ 
diation field (Fpuv) distribution is obtained via integration over the 
912 — 2066 A (6.0 — 13.6 eV) range at every grid point. Isotropic 
scattering is included in the model. In this way, thanks to the curved 
cavity rims, direct irradiation of the e nvelope by the proto star is 
also simulated (as illustrated in fig. 5 of lSruderer et alj2009l) . Each 
setup was computed 10 times and the median Ffuv was adopted 
for each grid cell, in order to minimize the numerical noise stem¬ 
ming from the Monte-Carlo nature of RADMC, while remaining 
computationally economic. 


2.2 Complex organic molecule chemistry 


The evolving composition of the system is computed with a two- 
phase (gaseous and solid) chemistry code designed to solve for 
the abundances of species in time. The reaction rates depend on 
the physical condit ions and the chemical network, as detailed in 
IWalsh et 'H] <201 4alf9). The ne twork has also been used previously 


Drozdovskava et ^ <2014 . with only minor updates incorpo¬ 


rated into the version used here. The chemical network is a com¬ 
pilation of the Rate 1 2 release of the UMI ST Database for As- 
trochemsitry (UDfA0, iMcElrov et alj 20131) and the Ohio State 
University (OSU) network Jcarrod. Weaver & Herbstll200^ . with 
a total of 668 species and 8764 reactions. 


^ http://www.ita.uni-heidelberg.de/~dullemond/software/RADMC/indi 
^ http://www.udfa.net 
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The chemical model includes gas-phase two-body asso¬ 
ciations, recombination of cations on grain surfaces, adsorp¬ 
tion, thermal desorption, CR-induced thermal desorption (i.e., 
spot heating), photodesorption, photodissociation and ioniza¬ 
tion in the gas and solid phases, grain-surface two-body reac- 
tions, and reactive desorption with an efficiency of 1 pe r cent 
lOarrod, Wakelam & Herbsll2007l : IVasvunin & Herbslll2013h . Two 

sources of FUV photons are accounted for: the protostar with the at¬ 
tenuated -Fpuv computed from RADMC, and the CRs, which pro¬ 
duce FUV photons through excitation of H 2 (the most abundant 
molecule). Within the code, Ffuv is converted to visual extinction 
(Av) via 

Av = t/3.02, (6) 

where r is the UV extinction dBohlin, Savage & Drakell978h given 
by: 


T = —In 


Fpuv 


TT X Bx in) dX X Rl/ (7?2 -h ^2) 


(7) 


where the denominator is the blackbody radiation for T* over the 
same FUV wavelength range as the numerator, with the inclusion 
of geometrical dilution, and where tt accounts for radiation stem¬ 
ming from one hemisphere t owards a point in the envelope. A sim - 
ilar concept was applied in IVisser, Doty & van Dishoeclj d201lh 
and IProzdovskava et ^ ll2014|). A CR ionization rate (Ccr) of 
5 X 10~^^ s~^ is used ( lDalgarn3l2006l) . which is slightly higher 
than the traditional dark cloud value ((^ism = 1.3 x s“^). 

FUV photons stemming from C R excitation of H 2 (with a flux of 
10 "* photons cm“^ s“^ for UsM- IShen et akllToOdh are produced in¬ 
ternally, and are thus available in the most dense and shielded re¬ 
gions. The elevated ^cr favors chemical complexity, although re¬ 
cent work suggests that the io nization rate may be as much as 3 or¬ 
ders of magnitude lower (e.g.,| Padovani. Ftennebel le & Gallil2013l : 
ICleeves. Adams & BerginlioisOcieeves et al.l2015h . It is assumed 
that CRs penetrate the entire system uniformly, since even the max¬ 
imum gas column density (~ 1 — 2 g cm“^) is too low for CR 
attenuation. Furthermore, magnetic fields are not modelled in this 
work; and stellar winds are expected to be importa nt only in the 
innermost regions dCleeves. Adams & Berginll2013l) . External ir¬ 
radiation is not considered, since it strongly varies from source 
to source and is expected to affect only the very outermost en¬ 
velope. FI 2 , CO and N 2 self- and mutual shielding is accounted 
for in the calculation of the respective photodis sociation rates 
dvisser. van Dishoeck & Blackl2009l ;lLi et al j2013h. The rn ain dif¬ 
ference between the chemical setup in IProzdovskava et alj ( I20l4) 
and the one here, is the us e of photorate! ^ for a 10 000 K rather 
than a 4000 K blackbody (| van Pishoec k, Jonkheid & van Flemerti 
bOOfitlvan Flemert & van Pishoeckiioosih 

To decrease the computational weight of the chemistry mod¬ 
ule, all grains are assumed to be 0.1 /rm in radius (a), which 
is near the peak of typical dust populations. The number density 
of grain surface sites (A^sites) is 1.5 x 10^® cm~^ with a barrier 
thickness of 1 A, under the assu mption of a rectangular barrier 
dUasegawa, Flerbst & Leundl 19921) . Grains are considered to be ei¬ 
ther neutral or negatively charged, with a total, temporally and 
spatially constant number density (a;grains) relative to nn. It is as¬ 
sumed that the gas temperatures equal the dust temperatures de¬ 
rived with RAPMC. This should only break down at the low¬ 
est densities ( 10 ® — 10 ^ cm~®), e.g., within the outflow cavities 


themselves (where the chemistry is disre garded in this wor k) and 
in a thin layer along the cavity walls 1 Visser et al.l 12013) . Fur¬ 
thermore, gas temperature is of limited influence on the abun¬ 
dance of complex organics, because grain-surface associations and 
thermal desorption both depend solely on the dust temperature. 
Grain-surface reaction rates are computed assuming the Langmuir- 
Hinshelwood mechanism only and using the rate-equation method 
lUasegawa, Herbst & Leun3l 19921) . 

Currently, in the chemical network, the formation of complex 
organic molecules relies primarily on associations of radicals upon 
FUV irradiation and the warm-up of icy mantles. At the lowest tem¬ 
peratures, hydrogenation is assisted by quantum tunneling of H and 
H 2 between grain surface sites, which becomes faster than classi¬ 
cal thermal hopping. Quantum tunneling through reaction barriers 
is also allowed for these two species. The mobility of any species 
on the grain surface depends on its binding energy (iSdes, also called 
its desorption energy). Here, the diffusion barrier (iSdiff) is set by: 


77diff — 0.3 X F/des, 


(8) 


as in iHasegawa. Herbst & Leu^ ( 1 19921) . This allows higher mobil¬ 
ity at low tem peratures (~ 15 — 20 K) in comparison to the factor 
of 0.5 used in iGarrod & Herbsj (l2006h . Off-lattice kinetic Monte 
Carlo techniques suggest that the value is species-dependent, e.g., 
0.31 for CO and 0.39 for C O 2 have been determined on t he sur¬ 
face of crystalline water ice iKarssemeiier & CupDen|[^014ll - val¬ 
ues which are in line with that assumed in this work. The set of 
adopted iJdes values is the one supplied with Rate 12, with the ex¬ 
ception of water, fo r which a higher va lue of 5773 K for pure water 
ice is used instead jFraser et al.ll200ll) . These values can vary sig¬ 
nificantly depending on the surface composition, which may vary 
with the formation sequence of icy mantles and eventual segrega¬ 
tion of solid species within the mantle. 

FUV photons are responsible for photodissociating molecules 
and producing radicals, which may go on to make larger species. 
Radicals can be formed directly in the ice or in the gas phase, 
followed by subsequent adsorption or gas-phase reactions. How¬ 
ever, FUV photolysis is a delicate balance. A strong FUV radia¬ 
tion field breaks molecules apart and favorably produces photofrag¬ 
ments (like CH 3 and C 2 H 5 ). In this work, an elevated ^cr favors 
radical production in strongly a ttenuated regions (Ay > 5 mag) by 
CR-induced FUV photons (e.g., [^rakines. Schutte & Ehrenfreundj 
ll996l : IPalumbo et^l2008l), which has been seen to e nhance com¬ 
plexity in experiments jMaitv. Kaiser & Jonesll2015l) . Elsewhere, 
the stellar FUV field dominates. Other mechanisms for produc¬ 
ing radicals have been postulated in the literature and supported 
by ex periments, e.g., low-energy electron radiolysis jBoamah et alj 
l2014h . but they are not taken into account in this work. 

The assumed grain properties are highly simplified in this 
chemical setup. Grain sizes and their distribution change with time, 
although theory su ggests that this is most prominent for proto¬ 
planetary discs (see iTesti et al.ll2014l for a review; however there 
is also accumulating obs ervational evidence for grain growth in 
the embedded phase , e.g., Jgrgensen^ ^l2009l; iRicci et al J|20 1 d; 


iLommen et al 1120 id; Ip agani et alJl20ld;lMiotello et a 1J|20 i 4) . The 

grains are also not necessarily spherical, and their surfaces may be 
highly irregular. All these parameters will change the surface area 
available for grain-surface chemistry and the number of binding 
sites that are available per unit volume. Under the current assump¬ 
tions, the key parameter is the total density of sites (nsues), given 
by: 


http://home.strw.leidenuniv.nl/~ewine/photo/ 


'^sites —- 47r X Cl X -A^sites X Xgrains X 7l }\. 


(9) 
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Figure 1. The physical sttucture of the system, from top to bottom: the gas 
density, nn (cm~®), the dust temperature, Tjust (K) and the visual extinc¬ 
tion, Ay (mag). The outflow cavity wall is shown with a black curve. Data 
points only within a fixed radius (16 100 AU) of the protostar are consid¬ 
ered. 


Here, Xgrains = 2.2 x 10“^^ is assumed. Decreasing this number by 
a factor of ~ 2, for example, results in abundances of gaseous and 
solid species changing by a factor of ~ 2 — 3 non-linearly. 

Finally, time is also an important parameter. Grain-surface 
radical-radical association reaction rates are highly temperature- 
dependent and may require a long time to build up appreciable 
product abundances. Thermal desorption, on the other hand, can 
be very rapid if the dust temperature is sufficiently high. Radicals 
may thus be lost to the gas phase before a reaction can take place. 
In this work, the chemical evolution of the system is computed up 
to 3 X 10® yr, which correspond s to the upper limits on the age of 
embedded sources (see table 1 in lPunham et alj[20l4) . 


3 RESULTS 

3.1 Physical structure 

The physical structure of the system is displayed in Fig.[T] In the top 
panel, the number density of H nuclei follows the adopted power 
law throughout the envelope, nu oc (equation[3]and Table[T](. 
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Figure 2. The abundance of methanol in the solid (left column) and gas 
(right column) phases at six different time steps across the envelope-cavity 
system. The outflow cavity wall is shown with a black curve. White cells 
correspond to either being outside of the area being considered or to having 
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Figure 3. Same as Fig. [2] but for formic acid. 


Figure 4. Same as Fig.[2| but for methyl formate. 
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Figure 5. Same as Fig. [2] but for acetaldehyde. 


Figure 6. Same as Fig. [2] but for dimethyl ether. 
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Table 2. Abundances at the end of the prestellar core phase”. 


Species 

Name 

n (Xgas) /m 

n (Xice) /riH 

CO 

carbon monoxide 

3.5 

X 

10 “® 

5.6 

X 10“® 

H 2 O 

water 

5.7 

X 

10 “® 

2.0 

1 

0 

X 

H 2 CO 

formaldehyde 

6.0 

X 

10 “® 

1.4 

X 10“® 

CH 2 CO 

ketene 

6.7 

X 

10 “® 

1.7 

00 

1 

0 

r-H 

X 

CH 3 OH 

methanol 

1.3 

X 

10 “®° 

9.2 

X 10“® 

HCOOH 

formic acid 

1.1 

X 

10 “®° 

4.4 

X 10“®° 

HCOOCH 3 

methyl formate 

3.2 

X 

10 “®® 

4.4 

X 

0 

1 

CA 

CH 3 CHO 

acetaldehyde 

2.1 

X 

10 “° 

1.1 

X 10“° 

CH 3 OCH 3 

dimethyl ether 

3.1 

X 

10 “®® 

3.2 

X 10“®° 

C 2 H 5 OH 

ethanol 

5.1 

X 

10 “®° 

4.3 

X 10“®° 

CH 3 COOH 

acetic acid 

2.7 

X 

10 “°® 

5.5 

X 10“®® 

HOCH 2 CHO 

glycolaldehyde 

8.6 

X 

10 “°® 

2.4 

X 10“®'^ 


“ The physical conditions of the prestellar core phase are assumed to be 
nii = 4 X IC* cm“®, Tjust = 10 K, Tfuv = 0 erg cm“^ and 
t = 3 X 10® yr. 


The outflow cavity has a fixed low density of 2 x 10^ cm“® (equa¬ 
tion |2l(- The ellipsoidal cavity wall, as specified by equation [T] for 
the parameters given in Table [T] separates the zone carved out by 
the outflow from the envelope. The middle panel shows the dust 
temperature as obtained with RADMC. The cavity is the warmest 
region with Tdust > 100 K within the inner ~ 300 AU of the pro¬ 
tostar, as a result of its low density and thus low extinction. The 
envelope is predominantly cold, Tdust < 25 K. Only a thin strip of 
material along the cavity wall (for 2 ; < 8000 AU) attains lukewarm 
conditions, ~ 30 — 40 K. Along the midplane, the Tdust = 100 K 
threshold is exceeded at i? « 38 AU. This, in combination with 
our exclusion of the cavity itself, implies that the traditional hot 
corino is not treated in this work. The bottom panel displays the 
visual extinction, which is computed with equations|^and|7]based 
on the stellar FUV radiation field calculated with RADMC. Again 
due to its low density, the cavity is subject to the strongest irradi¬ 
ation, Am < 2 mag. The envelope zones closest to the radial axis 
are fully shielded from stellar FUV photons, A\ > 10 mag, thanks 
to the large dust column density separating them from the proto¬ 
star. The rest of the envelope is subject to weak FUV irradiation. 
Am < 3 mag. 

The inner 35.9 AU around the protostar are excluded from 
this setup (Table [T), where the physical conditions need to be ad¬ 
justed for the presence of a protoplanetary disc. The RADMC cal¬ 
culations are performed on a much finer grid than that used for 
plotting in Fig. [T] which is the adopted grid for time-consuming 
chemical computations (4290 points in total). All cells are treated 
individually in the subsequent chemical calculations, i.e., as in a 
static model (in contrast w ith a dynamic evolutionary setup, as in 
IProzdovskava et ^l2014h . This maximizes the time spent in each 
temperature regime. 


3.2 Initial abundances 


The output of a single-point model of the prestellar core phase 
under constant physical conditions is assumed to be representa¬ 
tive of the envelope material and therefore, is adopted as the set 
of initial molecular abundances. This s chem e was also followed 
by Visser. Doty & van Dishoec^ ( 1201 ih and IProzdovskava et ^ 
(I2OI4I) . The assumed physical parameters of the prestellar core are 
nn = 4 X 10"* cm~®, Tdust = 10 K, neglig ible stellar FUV flux 
and an age of 3 x 10® vr. I Andre et al.l ( 1201 4 I) (in their section 3.3) 


estimated that the average lifetime of a starless core with a typi¬ 
cal density of ~ 10"^ cm“® is ~ 10® yr. A selection of the ob¬ 
tained abun dances is provided in Table Dark core observa tions 
of L1689B dBacmann et alJlToi2h . L1544 dVastel et al. 2014) and 
observa tions towards the ‘core’ position in Bl-b of Oberg et ^ 
d2010bll (and data from ICemicharo et al.l2012|) . derive column den¬ 
sities on the order of ~ 10^^ cm~^ for ketene, formic acid, methyl 
formate, acetaldehyde and dimethyl ether. If a typical H 2 column 
density of ~ 10^® cm“^ is used, then the observed abundances of 
these species in prestellar cores are ~ 10~^^ within an order of 
magnitude. A comparison with the values in Table reveals that 
from the aforementioned list only gas-phase methyl formate and 
dimethyl ether are underestimated in this prestellar core model. 
On the other hand, the modelled abundances of simpler ices, like 
water, carbon monoxide and methanol, are clos e to the observed 
ice abundances against background stars (e.g., Kngz_etaL 20051 : 


iBoogert et alJl201 ll. 1201 3l: iBoogert, Gerakines & Whittell2015l). 


3.3 Abundance maps 

Figs I2I6I show the abundance maps of five different molecules at 
six different time steps. The left six panels display the solid phase 
and the right six panels show the gas phase. For completeness, the 
abundance maps of five additional molecules (water, formaldehyde, 
ketene, ethanol, acetic acid) are provided in Appendix IaI 1 and gly- 
colaldehyde is discussed in Section [3.4.3t . 


3.3.1 Methanol 

Fig.lH shows that the methanol ice abundance remains predomi¬ 
nantly unchanged from its initial abundance (9.2 x 10' -®, TableO 
for the first 10® yr. With time, the abundance drops outwards from 
the cavity wall and eventually methanol ice is abundant (~ 10“®) 
only in the lower half of the envelope. This is the effect of photodis¬ 
sociation of methanol ice by stellar FUV photons (recall the Am 
map from Fig.[T). At 10"^ yr, the upper cavity wall layer displays 
a small enhancement (within an order of magnitude). Stellar FUV 
photons do not only destroy methanol ice, but they also produce 
CH3 and OH radicals by photodissociating other species (like H2O 
ice), providing an additional formation channel for methanol ice 
(besides sequential hydrogenation of CO, which is efficient across 
the entire envelope). This enhancement along the wall survives un¬ 
til 10® yr, when the more distant zone (at 2 ~ 12 000 AU and 
7? ~ 14 000 AU) starts to show a lower abundance (~ 10“^) by 
an order of magnitude, indicating favouritism towards destruction 
in that region. 

Another methanol ice-rich zone is formed around 2 ~ 500 AU 
and R ~ 5000 AU. It will be shown in subsequent maps (Fig. [Hi 
that this small area is rich in all complex organic ices, and for that 
reason this area is referred to as the ‘complex organic molecule 
torus’. Methanol, being a key parent species, is also abundant here. 
As it is converted to larger species, the enhancement is subdued 
at 10® yr. At longer time-scales, ~ 3 x 10® yr, CR-induced FUV 
photons start to dominate and photodissociate larger molecules to 
remake methanol, leading to the regaining of the enhancement and 
a resetting of the ice chemical complexity. 

The distribution of gas-phase methanol is closely related to 
the abundance of methanol ice. In the regions furthest away from 
the cavity, where methanol ice is abundant (~ 10“®) for all times, 
the abundance of methanol gas is low (< 5 x 10“^®). The en¬ 
tire envelope is colder than the thermal desorption temperature of 
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methanol 90 K), and there are no efficient gas-phase routes to 
methanol, therefore the amount of gaseous methanol relies purely 
on non-thermal desorption. In the furthest regions, extinction is 
high, so photodesorption with stellar FUV photons is insignificant. 
Gas-phase methanol is obtained through weak photodesorption by 
CR-induced FUV photons and reactive desorption stemming from 
hydrogenation reactions with CO. The cavity walls, on the other 
hand, light up in gaseous methanol for the first 10® yr. This 
corresponds to the zone showing the methanol ice enhancement. 
There, enough stellar FUV photons penetrate to allow more effi¬ 
cient photodesorption, but more importantly, to dissociate methanol 
ice. Upon the subsequently recombination of the photoproducts, 
reactive desorption is assumed to occur in 1 per cent of such reac¬ 
tions, enhancing the gas-phase methanol abundance. 

3.3.2 Formic acid 

Formic acid has an initial abundance of ~ 10“^° in both phases. 
There are several key grain-surface pathways leading to its forma¬ 
tion, specifically hydrogenation (H -l- COOK) and associations of 
OH either with H 2 CO or HCO. Solid formic acid is clearly en¬ 
hanced along the cavity wall at all time steps in Fig. [3] This occurs 
due to elevated radical-radical association rates thanks to stellar 
FUV photons, as was the case with methanol ice. The enhancement 
moves further away from the cavity rim with time with the separa¬ 
tion between the two being largest 4000 AU) at highest « values 
(~ 10 000 AU), which is a result of the formation-destruction bal¬ 
ance shifting. The zone with the highest abundances (> 10“®) is 
maximized at 10® yr. Gas-phase formic acid illuminates the cavity 
wall primarily at the earlier times 10® yr). Its abundance peaks 
(> 5 X 10“®) at 10^ yr, before that of the ice. This corresponds 
to the time when radical-radical chemistry is most active. On the 
longest time-scales, CR-induced FUV photons destroy formic acid 
in both phases. 

3.3.3 Methyl formate 

Fig. a shows the distribution of methyl formate. Its ice reaches 
abundances of ~ 10 “® along the cavity wall, which is a large in¬ 
crease from its initial abundance of 10“^®. Peak abundances 
> 10 “^ are reached at 10 ® yr in the complex organic molecule 
torus. Gaseous methyl formate exceeds the ~ 10“^^ abundance 
level at early times along the cavity wall, but mostly remains at a 
low abundance (< 10“^^). The enhancement of the ice in the com¬ 
plex organic molecule torus is not reflected in the gas-phase abun¬ 
dance, because at its location, there are no efficient routes to liber¬ 
ate the icy mantle species. This was also the case with methanol and 
formic acid in Figs[2and[3 respectively. The formation of methyl 
formate predominantly occurs either via the association of methoxy 
with HCO or with CO followed by hydrogenation. 

3.3.4 Acetaldehyde 

Acetaldehyde ice remains at its initial abundance of ~ 10“® 
throughout the envelope. With time, photodissociation by stel¬ 
lar FUV photons depletes it in the vicinity of the cavity wall 
(Fig. [3. An appreciable enhancement is built up in the torus with 
an abundance of 10“^ — 10“®. Gaseous acetaldehyde appears 
to be highly decoupled from its ice due to a strong enhancement 
(> 10 “®) along the wall for < 10 ® yr, which is not seen in the 
solid phase. This feature appears due to a single gas-phase reaction. 


namely O - 1 - C 2 H 5 (with a rate coefficient of 1.33 x 10“^® cm® s“® 
upon extrapolation from iTsang & Hampsoiil Il986h . Atomic O is 
highest in abundance towards zones with the stronger irradiation, 
i.e., towards the cavity, while C 2 H 5 is most abundant towards zones 
with the least irradiation, i.e., towards the radial axis. This results 
that in the ‘overlap region’, the cavity wall layer, the reaction is 
very efficient at producing acetaldehyde. Grain-surface production 
of acetaldehyde proceeds via CH 3 - 1 - HCO and CH 3 - 1 - CO followed 
by a hydrogenation. 


3.3.5 Dimethyl ether 

For the first 10® yr, dimethyl ether strongly highlights the cav¬ 
ity wall with solid and gaseous abundances higher than 10 “^ and 
10“®, respectively (Fig. n. This is the result of active grain-surface 
chemistry and concomitant reactive desorption. At longer times 
(> 10 ® yr), the enhancement is destroyed due to the domination of 
photodissociation by stellar FUV photons. However, at those times 
dimethyl ether ice is built up in the torus at an abundance of 10 “^, 
which is hidden from the gas phase for the above-stated reasons. 
The only grain-surface pathway towards dimethyl ether included in 
the network is the association of CH 3 with methoxy. 


3.4 Parameter study 

3.4.1 Reactive desorption 

During the above analysis of abundance maps, the importance of 
reactive desorption is readily seen, even with a mere 1 per cent 
efficiency. This has also been reported in previous publications. 
IVasvunin & Herbsd ( 120 13h saw variations of several orders of mag¬ 
nitude in gas-phase abundances and differences of 10 per cent in 
ice abundances for a dark core model (T = 10 K, nn = 10® cm“®. 
Ay = 10 mag, t = 10 ® — 10 ® yr ) upon varying r e active desorption 
efficiency from 0 to 10 per cent. IWakelam et alj ( l2014h noted that 
gas-phase water abundances can change by a factor of 10 upon ex¬ 
clusion of reactive desorption (for various single point models with 
T = 15-30K,nH = 2xl0'‘-2xl0® cm“®,Av = 2-4mag,f = 
10^ — 1 0® yr). Experiments con straining this process are sc arce. 
Recently. iDulieu et^ j2013l) and iMinissale & DulieJ ( l2014h sug¬ 
gested that the efficiency of reactive desorption can vary greatly 
per reaction and per surface, from very efficient on bare grains to 
completely inefficient as soon as the first monolayer is built up. 

Fig. El shows the ratio of the methanol abundance with reac¬ 
tive desorption switched off to that with an efficiency of 1 per cent, 
i.e., comparing with Fig.|3 Gas-phase abundances decrease over 
almost the entire envelope for all times once reactive desorption is 
excluded. The decrease is 1 — 3 orders of magnitude for the largest 
portions of the envelope, and is explained by the exclusion of the 
most efficient mechanism for the population of the gas phase by the 
methanol-rich ices. In strongly irradiated areas. Ay < 3 mag, the 
abundance of solid methanol is as much as 3 orders of magnitude 
higher when reactive desorption is not included, which is an ice de¬ 
struction mechanism. However, in zones where CR-induced FUV 
photons dominate, Ay > 3 mag, there is less solid methanol when 
reactive desorption is excluded. Reactive desorption elevates the 
gas-phase abundance not only of methanol and complex organic 
molecules, but also of simpler species, like CH 4 . Gas-phase ion- 
molecule reactions then lead to an enhancement of gas-phase radi¬ 
cals, most notably CH 3 . Higher gaseous radical abundances imply 
greater availability of these species on the grains (transiently) as 
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Figure 7. The ratio of the methanol abundance with reactive desorption 
switched off to that with an efficiency of 1 per cent. Three time steps are 
shown with the solid-phase and the gas-phase ratios in the left and right 
columns, respectively. The outflow cavity wall is shown with a black curve. 
A certain colour corresponds to a range of values between the two labeled 
bounds, e.g., yellow corresponds to any value between 1 and 10. 


Figure 8. The ratio of the methanol abundance with ice photodissociation 
(by stellar and CR-induced FUV photons) switched off to that with it on. 
Three time steps are shown with the solid-phase and the gas-phase ratios in 
the left and right columns, respectively. The outflow cavity wall is shown 
with a black curve. 


well, thus allowing more efficient formation of methanol via the 
OH + CHa route. 

Acetaldehyde and dimethyl ether follow a trend very similar to 
that seen with methanol in Fig.|7] Formic acid and methyl formate 
differ by not showing a decrease in the zone with high extinction 
when reactive desorption is switched off. This is linked to them 
being less reliant on the availability of CH 3 . 


3.4.2 Ice photodissociation 

Another poorly constrained process is the photodissociation of 
solid species. UV photons can penetrate as deep as 100 mono- 
layers into the icy mantle; however, what happens with the 
photofragments thereafter remains unclear. They may recom¬ 
bine immediately or diffuse away through the ice, if they 
have sufficient translational energ y after photodissociation (e.g., 
lAndersson & van Dishoeckl 1^081) . In the current setup, it is as¬ 
sumed that the entire mantle may be dissociated; however, by mak¬ 
ing all those radicals equally available for further grain-surface re¬ 
actions, the chemistry of the solid phase is potentially overesti¬ 


mated, because the diffusion rates within the bulk ice mantle are 
likely slower than those across the surface. 

To test the significance of ice photodissociation by stellar and 
CR-induced FUV photons, a control simulation is carried out with 
those processes turned off. Photodesorption and reactive desorption 
are still included. Fig.[ 8 ]shows the ratio of the methanol abundance 
with ice photodissociation excluded to that with it included, i.e., 
again comparing with Fig. The amount of methanol ice is in¬ 
creased roughly by an order of magnitude in zones with the highest 
extinction {Am > 3), and in the cavity wall by 3 orders of magni¬ 
tude. This corresponds to where photodissociation is dominated by 
stellar UV photons and is the primary ice destruction pathway. For 
t ^ 10 ® yr, the gas-phase methanol abundance is reduced, since 
less ice photodissociation implies less radical recombinations, and 
thus less gas-phase methanol due to reactive desorption. However, 
at later time steps and in the cavity walls, a much larger abun¬ 
dance of ice leads to more gaseous methanol via photodesorption, 
CR-induced thermal desorption, and reactive desorption associated 
with slow formation pathways. At the final time step, the gas-phase 
methanol abundance starts to decrease in the most shielded zone, 
since CR-induced FUV photons are still photodissociating in the 
gas-phase, but efficient mechanisms (in particular, reactive desorp- 











































The chemistry of outflow cavity walls exposed 11 


tion, in comparison to the fiducial setup) to replenish the gas from 
the ice are absent. 

Other complex organics, i.e., formic acid, methyl formate, ac¬ 
etaldehyde and dimethyl ether, show a decrease by several orders 
of magnitude in both phases once ice photodissociation is switched 
off. In this control simulation, the radical production is limited to 
the gas phase, while key radical sources, such as methanol, are most 
abundant in the solid phase. The reduction of radicals drastically 
obstructs efficient formation of complex organic molecules. Only 
along the cavity wall, are the abundances higher than fiducially. 
Ice photodissociation dominates there and in the standard case ef¬ 
ficiently destroys complex organics. 


points in the COM torus (along 2 = 500 AU and for 3756 < 
R < 6006 AU) gives a ratio in the ice of 0.84 fiducially and of 
0.04 after the route via glyoxal is introduced. The models indi¬ 
cate that there is more glycolaldehyde than methyl formate, while 
observations suggest the contrary. This may be explained by the 
fact that the observations probe predominantly the hot core re¬ 
gions , where gas-phase reactions become important (iTaguet et alj 
UOlSh . They may drive chemistry that is different from that oc- 
curing in the sol id phase. This ma y include gas-phase reactions 
with formic acid dTaguet et aklbOlsT) . production of methoxy from 
OH-t-CHaOH and hydrogen-abstract ion reactions with F and Cl 
ISalucani, Ceccarelli & Taouetll201^ . 


3.4.3 An additional route to glycolaldehyde 

The chemistry explored in this work is limited by the reactions and 
the species included in our chemical network. In order to inves¬ 
tigate this limitation, an additional formation pathway for glyco¬ 
laldehyde is introduced, namely: 

HCOice -b HCOice ^ OCHCHOice, 


„ HOCHCHOice „ 

OCHCHOice 4 , 4 HOCHaCHOice. 

OCHaCHOice 

This scheme was sugg ested bvlWoods et alJ 2013h and experimen¬ 
tally proven for ices bv iFedoseev et al | |2015 ). In total, 68 reactions 
are added to the existing chemical network and are detailed upon 
in Appendix IB] 

The same initial conditions are used as in the fiducial case 
(which also holds for all other test cases discussed), so that the de¬ 
pendence of the results purely on the inclusion of the additional 
pathway is explored. The abundance of gaseous and solid glyco- 
ladehyde increases by an order of magnitude across the system 
when the additional formation route is included, as Fig. shows 
(the difference is larger in areas where the abundance of glyco¬ 
laldehyde is very low (n(HOCH 2 CHOice)/n-u < 10 “^^), but that 
is not significan t ). Thi s is consistent with the upper limit derived 
by IWoods et alj (l2013h . where destruction of glycolaldehyde and 
competitive routes are not considered. Adding an extra route for 
forming a molecule is expected to boost its abundance. The boost 
is larger in highly extincted and, in turn, cold (Tdust < 20 K) zones, 
where hydrogenation dominates over radical-radical associations. 
This enhancement also prolongs the lifetime of glycolaldehyde, and 
it is still found at appreciable ice abundances (~ 10~®) in the COM 
torus at 10® yr, which was not the case previously. All other species 
show differences by a factor of a few or less (any modification of 
a chemical network results in non-linear effects on all species due 
to competing pathways, e.g., the inclusion of HCO - 1 - HCO com¬ 
petes with HC O + H, HCO + CH 3 O, e tc.), which is consistent with 
the findings of iFedoseev et alj l l2015h using Monte Carlo models. 
This exercise shows that there may remain important routes to the 
formation of complex organics that have not yet been included in 
models. 

Glycolaldehyde and methyl formate (and acetic acid) are iso¬ 
mers and may be intriguing probes of the physics and chemistry 
of a system. Recent observations have estimated the methyl for- 
mate to glycoladehyd e ratio to be ~ 13 for IRAS 16293-2422 
jjdrgensen et alj 2of^ Coutens et al.l 2015h . 12 — 20 for NGC 

1333-IRAS2A dCoutens et al.ll2015l:lTaouet et al.ll2015l) and - 10 
for NGC 1333-IRAS4A i Taouet et all 20151) . Averaging over 10 


3.4.4 Stellar luminosity 

Young Class 0 and I protostars cover a range of luminosities 
from fractions of L,^ to more than an order of magnitude higher 
llEvans et alJl2009l : lDunham et al.ll2014l) . The dependence of these 
models on the luminosity is investigated by computing the chem¬ 
istry for the case of 1 Lq and 15 Lq. Our adopted luminosity of 
35.7 Lq for NGC 1333-IRAS2A is on the higher side of the bulk of 
the low-mass protostars. Since the density is fixed, and both Ffuv 
and the geometric dilution of the blackbody radiation scale with L* 
(see equation|7]l, the A\ map remains unchanged for different lumi¬ 
nosities. However, the number of FUV photons reaching a certain 
grid cell does change. Fig. lClI in Appendix O shows the difference 
in the dust temperature between the fiducial 35.7 Lq model and 
the two with lower luminosities. In the case of 15 Lq, the envelope 
dust temperatures are cooler by at most 5 K. For 1 Lq, the differ¬ 
ences are larger, and the cavity wall can be as much as 10 — 15 K 
colder than in the original case. 

Fig. IC2I shows the methanol abundance in the case of 15 Lq 
and 1 Lq in comparison to the fiducial setup. For solid and gaseous 
methanol (with both phases displaying strong coupling), the abun¬ 
dances for 15 Lq and the fiducial run are very similar with dif¬ 
ferences predominantly within 1 order of magnitude. Only small 
regions show an increase in methanol by 2 — 3 orders of magni¬ 
tude in the 15 Lq case. The enhancement along the cavity wall is 
explained by the shift of the temperature regime most favorable for 
grain-surface chemistry due to enhanced mobility prior to thermal 
desorption of radicals and the decrease of FUV photons making 
radical production inefficient. Other species (formic acid, methyl 
formate, acetaldehyde and dimethyl ether) predominantly show a 
decrease in abundance by an order of magnitude, leading to the in¬ 
crease in simpler species, like methanol and formaldehyde. 

For the 1 Lq case (left column of Fig. IC2L the trends seen 
with 15 Lq are dramatized. The formation of complex organ¬ 
ics is significantly impeded with gaseous and solid abundances of 
~ or lower in the entire envelope. This results in an increase 

of solid, and in turn gaseous, methanol, by as much as 3 orders of 
magnitude in a broad zone along the cavity. In conclusion, the orig¬ 
inal model with the highest luminosity of 35.7 Lq has been identi¬ 
fied as most efficient for chemical complexity within this paradigm, 
while lower luminosities support stronger enhancement or retain- 
ment of simpler species, such as methanol. The morphology of the 
region, including the highlighting of the cavity wall, is preserved 
upon varying L*, although the sizes of these regions change as do 
the peak abundances. 
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Figure 9. The abundance of glycolaldehyde in the solid (upper six panels) and gas (lower six panels) phases at three different time steps across the envelope- 
cavity system. The right column corresponds to the fiducial model. The left column is the model output upon the inclusion of the additional formation pathway 
via glyoxal. The outflow cavity wall is shown with a black curve. White cells con'espond to either being outside of the area being considered or to having 
values outside of the range of the colour bar. The range of the gas colour bar is different from the range of that of the ice. 




















































The chemistry of outflow cavity walls exposed 13 


Table 3. Outflow cavity parameters (additional model setups). 


Parameter 

Units 

Narrow 

Wide 

(2cav 

AU 

4.5 X 10"* 

3.4 X 10"^ 

bcav 

AU 

2.0 X 10® 

8.0 X 10® 

a{z = lOOOAU) 

deg 

45 

125 

a (z = 10 OOOAU) 

deg 

14 

59 


3.4.5 Outflow cavity full opening angle 

The shape of the outflow cavities and the cavity full opening angle 
set the flux of stellar photons penetrating the envelope, as men¬ 
tioned in Section [231 To test the importance of the full opening 
angle, two additional simulations have been executed, one with a 
narrower and another with a wider full opening angle (Table O. 
Fig. lC3l in Appendixlclshows the difference in the dust temperature 
between the fiducial case and the two additional models, as well as 
the ratios in extinction. The temperature variations are largest in the 
zone that in one model is part of the cavity and not in the other, as 
well as in the cavity wall. The variations can exceed 20 K close to 
the star, and are a result of the change in the density from the ra¬ 
dial profile of the envelope to a fixed low value of the cavity. When 
comparing the three extinction maps, the pattern from narrow to 
fiducial to wide angle is that of an opening flower blossom. In the 
ratio plot, that is reflected by a large scale decrease in extinction 
of 7 — 25 times in the cavity wall when going from the small cav¬ 
ity to the fiducial setup. When switching from the fiducial to the 
large cavity, a further decrease by a factor of 7 — 10 is seen in an 
area immediately after the zone with the previous weakening. Fur¬ 
thermore, there is also a change in extinction in the zone that is 
switching from envelope to cavity area, as was the case with the 
dust temperature. 

Fig. IC4I shows the methanol abundance in the case of a wider 
and a narrower full opening angle in comparison to the fiducial cav¬ 
ity. The most prominent change is the enhancement by 2 — 3 orders 
of magnitude ~ 5000 AU away from the cavity wall, which is pre¬ 
ceded by a decrease of 2 — 5 orders of magnitude. This is due to the 
shift of the zone with efficient reactive desorption with the widen¬ 
ing of the cavity. Other regions and other species predominantly 
show variations within an order of magnitude. The morphology of 
the system and highlighting of the cavity walls is preserved. Only 
the COM torus, still rich in complex organic species, is slightly 
reduced in size. 

In the case of a small cavity (left column of Fig. lC4b . complex 
organic molecules are only produced in a narrow strip in the im¬ 
mediate vicinity of the cavity wall, because in all other regions of 
the system, very few FUV photons penetrate making radical pro¬ 
duction inefficient. Methanol benefits from this and an enhance¬ 
ment in both solid and gas-phase abundances is seen for a much 
thicker layer along the cavity wall. In the fiducial setup, stellar FUV 
photons would actively photodissociate methanol in that zone. In 
essence, the full opening angle of a cavity determines primarily the 
amount of FUV radiation entering the envelope. The morphologi¬ 
cal structure of the system remains comparable, but enhancements 
in methanol and other complex organic molecules shift angularly 
with the migration of the regions with optimal temperature and ex¬ 
tinction combinations. 



Figure 10. An illustration summarizing the key zones of the envelope- 
cavity system, and the major species and their phases therein. The motion 
of the cavity wall layer with time is indicated with two arrows. COM stands 
for complex organic molecules. 




Figure 11. Lifetimes of complex organic ices in the COM torus (at z = 
500 AU). The green bars indicate when the abundance of the specified 
molecule is above 5 X 10“®. The filled pink stars mark the time when the 
zone (in R) above this cut off value is maximized. The larger orange-yellow 
stars mark when the peak abundance is highest. 


4 DISCUSSION 
4.1 Morphology 

The abundance maps in Figs I2I6I suggest an evolving chemical 
morphology of the envelope-cavity system. The representation in 
Fig. QO] summarises all the information from the individual species 
in one global picture. First, there is the outflow cavity wall layer, 
which lights up in various species at different times. This re¬ 
gion corresponds to where photodissociation by stellar UV photons 
leads preferentially to the formation of various large solid species 
via radical-radical associations on the grain. Reactive desorption 
stemming from many recombinations of photoproducts leads to an 
enhancement of methanol and complex organics in the gas phase. 
Acetaldehyde is an exception to this, due to an efficient gas-phase 
pathway (O -l- C 2 H 5 ) leading to a large gas-phase abundance, but 
not in the solid p hase. Such exception al behaviour was also seen in 
and suggested bv ICodella et alj ( ItOISI) . 
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Secondly, there is the COM torus with the largest diversity 
of complex organic ices at abundances several orders of magni¬ 
tude higher than initially injected into the system. The largest, most 
complex species (e.g., dimethyl ether in contrast to formic acid) ap¬ 
pear later in time in the torus than in the cavity wall layer, but also 
survive much longer due to the mild radiation. This zone probes the 
time-scales of cool (15 — 25 K) quiescent (^v ~ 10 mag) grain- 
surface chemistry with CR-induced FUV photons being the only 
source of radicals. In essence, the cavity wall layer is the scaled- 
up version of the torus, based on the scaling of the FUV radiation 
fields. A stronger stellar FUV radiation field builds complex organ¬ 
ics faster and in a larger zone, but for a short period of prosperity. 
A weaker CR-induced FUV radiation field builds complexity more 
slowly, in a smaller zone that is shielded from powerful dissociat¬ 
ing stellar radiation. Both require cool dust temperatures (the dis¬ 
tinguishing feature from the prestellar phase), and both eventually 
destroy complexity via photodissociation. 

Finally, the rest of the envelope is where only relatively simple 
ices and gases reside. This pertains to species like water, methanol, 
formaldehyde, and ketene among others. Few changes happen in 
this zone compared with the initial conditions, because the temper¬ 
ature is low (< 15 K) and shielding is high {Aw > 3 mag). The 
chemical composition of the outer region remains similar to the 
prestellar phase. 

Several important lessons can be learned from the abundance 
maps, besides the chemical morphology of the system. As is seen 
for the COM torus, gases do not always reflect the composition of 
the ices. Effective non-thermal desorption mechanisms are needed 
to couple the two phases in the predominantly cold (< 100 K) enve¬ 
lope. Furthermore, there may be exceptional species, like acetalde¬ 
hyde, that have an efficient gas-phase route at particular physical 
conditions (like that of the cavity wall layer), which would result 
in an inferred ice abundance that is several orders of magnitude 
too high. Observations of solids are needed to accurately constrain 
the icy content, which remains a challenging task for solid com¬ 
plex organic molecules. Some species with isolated features in the 
5 — 10 /rm range can be searched for with future missions like the 
James Webb Space Telescope (JWST). 

4.2 Comparison with observations 

From Fig. [2 it can be seen that gaseous methanol illuminates the 
cavity wall as early as 10® yr with abundances of ~ 10“®. Formic 
acid, acetaldehyde, dimethyl ether and ethanol reach gas-phase 
abundances of~ 10“®, 10~®, 10~® and 10“®, respectively, in the 
cavity wall layer. These va lues are in rough agreement with the 
column densities derived bv lArce et alj l l2008h for formic acid and 
ethanol at the B1 position (-^ 10^® —10^“^ cm“^ giving abundances 
~ 10-10 _ 10“® cm~®), as well as those for acetalde hyde and 
dimethyl ether (~ 10^® cm“^) from lOberg et alj ( l201lh towards 
Serpens. This implies that thermal desorption (for Tdust ^ 100 K) 
and sputtering with the passage of shocks along cavity walls do not 
need to be invoked when interpreting observations on these scales 
(several thousand AU). These processes are more likely to be ef¬ 
ficient on smaller (several hundred AU) scales and shorter time- 
scales (~ 10® yr). The enhancements seen in the cavity wall layer 
rely heavily on reactive desorption, which dominates over much 
less efficient direct photodesorption. This highlights the need for 
further laboratory studies quantifying the efficacy of reactive des¬ 
orption as a mechanism for releasing complex organic ices into the 
gas phase at low temperatures. If it is inefficient, then shocks may 
once again be necessary to explain observations. Methyl formate 


does not fol l ow the observed abundances of lArce et al.l l l2008h and 
lOberg et"^ ( l201lh in our models, and in fact, is only efficiently 
made in the COM torus. 

By design of the chemical network used in this work, 
methanol plays a central role in the synthesis of larger, more com¬ 
plex species. However, it appears not only as a precursor, but also as 
a descendent (or postcursor), once CR-induced FUV photons dom¬ 
inate and destroy complexity at longer time-scales (> 3 x 10® yr). 
This scenario may be challenged as more atom addition reactions 
are included into the chemical network. 

Our models predict the COM torus - a solid-state sweet spot, 
which has not been seen observationally. The abundances in this 
zone may still change due to influences of the protoplanetary disk 
and other dynamical effects (Section l44t . For NGC 7538 IRS9 - 
a massive YSO, lOberg et'H] (1201 3h reported a change in chem¬ 
istry within 8000 AU of the protostar, which potentially links 
with efficient grain-surface chemistry. It remains unclear whether 
a high-mass source would be associated with more efficient path¬ 
ways towards enriching the gas phase from the ices rich in com¬ 
plex organics, in comparison to the low-mass case modelled in this 
work. Although, rece nt ALMA obse r vation s towards the high-mass 
IRAS 16547-4247 by iHiguchi et al.l j2015h suggest that methanol 
emission follows an hourglass morphology. 


4.3 Lifetimes of complex organics 

From the abundance maps it was seen that species peak at different 
times. Motivated by this, the ordering of complex organic species 
in time is shown in Fig. as they appear in the COM torus (a 
slice at 2 = 500 AU is assumed to be representative, which covers 
R G [1006, 16 100] AU). The green bars show when the abun¬ 
dance of an ice is above 5 x 10~®, i.e., when the solid phase is 
enhanced. The time of maximal spatial extent of this enhancement 
along R is marked with a filled pink star. When the peak ice abun¬ 
dance is reached, this is indicated with a larger orange-yellow star. 
The two do not necessarily overlap for all species considered. As 
was hinted earlier. Figs ID [S] H] and [TT] do not support the parti¬ 
tion of methyl formate with acetaldehyde into the group of ‘cold’ 
organics, contrasting dimethyl ether as a ‘hot’ species. Instead, ac¬ 
etaldehyde appears to have more in common with the morphol¬ 
ogy and lifetime of dimethyl ether. Yet, all three are formed in 
the COM torus at cool temperatures. A different organization of 
complex organic molecules by their lifetimes and times of peak 
abundance is attempted in Fig. This would imply that species 
like formaldehyde, ketene, methanol, formic acid, methyl formate, 
acetic acid and glycolaldehyde are ‘early’ types; while acetalde¬ 
hyde, dimethyl ether and ethanol are ‘late’ types. Fig.ll llindicates 
a chemical pattern that is consistent with the chemistry invoked in 
our model. Ketene, dimethyl ether and ethanol reach abundances 
> 5 X 10~® only at 10® yr, while acetic acid abundances do not 
stay high outside of the 10® — 10® yr range. This implies that sim¬ 
ple ices, like methanol and ketene, are converted into more com¬ 
plex organic i ces with the evolut ionary stage of the source (Class 
0 to Class I). lOberg et"^ l l2013h proposed that ketene, methanol 
and acetaldehyde belong to ‘zeroth-generation’ molecules, which 
is supported in this work for the former two species, but not for the 
latter (although the abundance of acetaldehyde peaks earlier than 
the maximization of the spatial extent of its enhancement). 
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4.4 Episodic accretion and dynamics 

A potenti al caveat to the modelli ng p resented thus far is episodic 
accret ion. Ivisser & BerginI ll2012h and ivisser. Bergin & Jdrgenserj 
120 ISh performed the first studies towards understanding the chem¬ 
ical implications of an accretion burst. The authors estimated that 
the effects should remain observable up to 10^ yr after such an 
event and that there should be an excess of gaseous CO. A factor 
10 increase in stellar luminosity (typical for a moderately strong 
burst) is expected to push the 100 K radius in the envelope out 
from ~ 38 to ~ 114 AU and the 20 K radius from ~ 3300 to 
~ 10 000 AU for our setup IVisser & BerginI[20121) . Thus strong 
heating of the dust leading to thermal desorption of all species, in 
turn transferring the complexity from the solid phase into the gas, 
only occurs for the inner envelope and does not influence the large 
scales studied here. However, if episodic accretion is accompanied 
by sufficiently strong FUV irradiation (stronger than the typical in¬ 
terstellar UV field), which would photodissociate complex organic 
molecules, complexity may be reset for regions of the system. Al¬ 
ternatively, if FUV irradiation and/or dust heating are mild, then 
complexity may be enhanced, especially now that the 20 K ra¬ 
dius lies much further out. The effects of episodic accretion must 
strongly depend on the exact physical properties of such an event, 
which is beyond the scope of this paper. 

In addition, the static model used in this work neglects the dy¬ 
namical effects of the system. On smaller scales, such as the inner 
500 AU studied in IProzdovskava et alj (l2014ll . dynamics defi¬ 
nitely play a crucial role due to the proximity to the growing pro¬ 
tostar and the expanding protoplanetary disc, which is also where 
material moves through regions with high temperatures. However, 
on much larger scales of several thousand AU that are focused on 
in this work, a static model may be sufficient as the physical con¬ 
ditions have a shal lower gradient with distance through the e nve¬ 
lope. Furthermor e,! Jorgensen, Schoier & van DishoecH l l2002tl and 
iKristensen et alJ |2012|) suggest that although Class 0 envelopes are 
infalling, by the Class I phase some envelopes are already showing 
expansion. 


5 CONCLUSIONS 

The envelopes encompassing young low-mass protostars are suc¬ 
cessively eroded by bipolar outflows, resulting in cavities that grow 
with time. These outflow cavities allow additional photons to pen¬ 
etrate the envelope, thus enhancing the heating of the dust and 
the irradiation of the material ther ein. Various observations (e.g., 
lArce et alll 20081 : lOberg et al.ll201 ih have hinted that this results in 
unique morphological features and gas-phase detections of com¬ 
plex organic molecules. This paper describes a 2D physicochemical 
model of such an envelope-cavity system (excluding a traditional 
hot corino), including wavelength-dependent radiative transfer cal¬ 
culations and a comprehensive gas-grain chemical network. The 
main results of our modelling are; 

(i) An envelope-cavity system has three distinct regions. Firstly, 
the cavity wall layer, which displays time-dependent enhancements 
in solid and gaseous abundances of methanol, formic acid, ac¬ 
etaldehyde, dimethyl ether and ethanol. Secondly, closer to the star, 
a ‘COM torus’ - a zone rich in complex organic ices. Thirdly, 
the remaining outer envelope, which is predominantly comprised 
of simpler ices and gases, like water, methanol, formaldehyde and 
ketene. 

(ii) Gases do not always reflect the ice composition - the COM 


torus is poor in gaseous complex organics, while having the high¬ 
est solid abundances of the entire system. The division between the 
two phases depends on the relative efficiencies of non-thermal des¬ 
orption mechanisms in different regions. 

(iii) Complex organic ices peak in the COM torus at different 
times suggesting unique molecular lifetimes. 

(iv) The gas-phase enhancements along cavity walls in select 
molecules are consistent with observations. In strongly irradiated 
regions (A\j < 3 mag), such as the cavity wall layer, photodisso¬ 
ciation in the solid phase is frequent. Subsequent recombination of 
the photoproducts leads to frequent reactive desorption. Although 
the poorly constrained efficiency is assumed to be a mere 1 per 
cent, the high frequency of reactive desorption leads to gas-phase 
enhancements of several orders of magnitude. Direct photodesorp¬ 
tion (using the laboratory constrained yield of 10“® molecules per 
photon) is found to be inefficient in comparison, indicating the need 
for experimental quantification of the rates for reactive desorption. 

(v) Photodissociation directly in the solid phase produces radi¬ 
cals in the icy mantle, the absence of which impedes the formation 
of complex organic molecules. For efficiency, a sufficient number 
of FUV photons needs to penetrate the envelope and elevated cool 
dust temperatures enable radical mobility on the grain surface. 

(vi) Consequently, a high stellar luminosity (~ 35 Lq) favours 
chemical complexity. A low stellar luminosity (-^ 1 Lq) results in 
negligible complex organic abundances (< 10“^^). 

(vii) Also, a sufficiently wide cavity (for example, 
a(z = 10 OOOAU) = 45°) is needed to directly irradiate the 
envelope. 

This work has shown that low-mass protostars have gas-phase 
enhancements in complex organic molecules when outflow cavities 
are present, as well as hidden zones that are rich in complex or¬ 
ganic ices. Future research will explore the intriguing extrapolation 
of this to the high-mass case, especially in light of recent ALMA 
observations of methanol e mission with an hou rglass morphology 
towards IRAS 16547-4247 jHiguchi et alj|2015h . 
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APPENDIX A: ADDITIONAL ABUNDANCE MAPS 

In addition to the five sets of abundance maps shown in the main 
text of the paper, another five molecules are presented here for com¬ 
pleteness. 
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Figure Al. Same as Fig.[2| but for water. 


Figure A2. Same as Fig. [2] but for formaldehyde. 
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Figure A3. Same as Fig. [2] but for ketene. 


Figure A4. Same as Fig.[2| but for ethanol. 
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APPENDIX B: DETAILS OF THE ROUTE TO 
GLYCOLALDEHYDE VIA GLYOXAL 

The details of the added route to glycolaldehyde via glyoxal are de¬ 
scribed in this appendix. In this work, the r oute is only added for 
the solid phase (unlike in IWoods et alJl^oT^ . since three-hody gas 
phase reactions require a radiative association mechanism, which 
is slow by nature. Associations of radicals are assumed to be barri¬ 
erless, which includes the reaction between the two formyl radicals 

(HCO) and the second hydrogenation. Hydrogenation of glyoxal 
(OCHCHO) has an adopted activation ene rgy (Ea) of 1100 K, 
based on quantum chemical calculations of lWoods et al.l ( l^l3ll . 

The chemical network used is thereby expanded by 10 
new species: gaseous and solid glyoxal, the cation of glyoxal 

(OCHCHO*^), the protonated glyoxal (OCHCH 2 O'''), the gaseous 

• • 

and solid two intermediate species (HOCHCHO, OCH 2 CHO), 
a joint cation (C 202 H;^), and a joint protonated form of the in¬ 
termediates (C202H^). The gas- phase chemistry of glyoxal is ex - 
tracted from the OSU network ( IGarrod. Weaver & Herbslll2008h . 
The gas-phase chemistry of the two intermediate species is as¬ 
sumed to be similar to that of other large hydrogenated radicals 
and glyoxal itself, from which the gas-phase reaction rate coeffi¬ 
cients are adopted. The dissociative reaction channels with He'*' 
and photodissociation channels are assumed to be 2HCO - 1 - H for 

HOCHCHO and H 2 CO -t HCO for OCHaCHO. Grain-surface 
chemistry of glyoxal and the two intermediates is not expanded 
further than the scheme provided above, with the exception of one 
competing reaction, namely: 

OCHCHOice 4 HaCOice -f HCOice, 

which is included in the OSU network (again Ea = 1108 K). 
The binding e nergy of glyoxal is taken to be 3200 K 
(OSU network, iGarrod. Weaver & Herbsj 1200 Sh and for the 
two intermediates calculated according to the prescription in 

IGarrod. Weaver & HerbstI ( 1200 Sh to be i?des(HOCHCHOice) = 

6080 K and i?des(OCH 2 CHO) = 3800 K (which accounts for 
the former molecule being able to H-bond due to the presence of 
an OH functional group). 


APPENDIX C: SUPPORTING PARAMETER STUDY 
FIGURES 


Figure AS. Same as Fig. [2] but for acetic acid. 
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Figure Cl. The difference in the dust temperature between the fiducial, 35.7 Lq model, and the two with smaller luminosities (15 Lq and 1 Lq). The 
outflow cavity wall is shown with a black curve. 


This paper has been typeset from a TgX/ KTgX file prepared by the 
author. 









22 Maria N. Drozdovskaya et al. 


n(CH30H)/nH 

15Le 


35.7Lm 


1.0x10^ 





1.0x10' 


R (10" AU) 


Figure C2. The abundance of methanol in the solid (upper nine figures) and gas (lower nine figures) phases at three different time steps across the envelope- 
cavity system. The left, middle and right columns correspond to the cases of 1 Lq, 15 Lq and 35.7 Lq (the fiducial model), respectively. The outflow cavity 
wall is shown with a black curve. White cells coiTespond to either being outside of the area being considered or to having values outside of the range of the 
colour bai‘. 
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Figure C3. The difference in the dust temperature between the fiducial model and the two additional models, one with a larger cavity (upper left panel) and 
another with a smaller cavity (upper right panel), as well as the ratios in extinction (lower left and right panels, respectively). The different outflow cavity walls 
are shown and labeled. 
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Figure C4. The abundance of methanol in the solid (upper nine figures) and gas (lower nine figures) phases at three different time steps across the envelope- 
cavity system. The left, middle and right columns correspond to the cases of a small, fiducial and large cavities, respectively. The outflow cavity wall is shown 
with a black curve. White cells coiTespond to either being outside of the area being considered or to having values outside of the range of the colour bar. 





























































